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Abstract 

Numerical differential equations studies the methods for 
finding numerical approximations to the solutions of 
differential equations, occur in many scientific disciplines, 
for instance in physics, chemistry, biology, and economics. 
Because of its various applications, this is often viewed as 
a discipline in and of itself. In this paper we develop a 
numerical simulator for solving ordinary differential 
equations (ODEs). This simulator is incorporated with a 
combination of Euler, modified Euler and Runge-Kutta 
second and fourth order methods. 

1. Introduction 

Very often the differential equations appearing in physical 
problems cannot be solved analytically. Thus it becomes 
imperative to discuss their solutions by numerical methods. 
In numerical method, we do not proceed in the hope of 
finding a relation between variables but we find the 
numerical values of the dependent variables for certain 
values of independent variables . 

The process of finding a derivative is called differentiation. 
The modern development of differentiation is usually 
credited to Isaac Newton (1643-1727) and Gottfried Leibniz 
(1646-1716), who provided independent and unified 
approaches to differentiation and derivatives 1 " 1 . In numerical 
analysis, numerical differentiation describes algorithms for 
estimating the derivative of a mathematical function or 
function subroutine using values of the function and perhaps 
other knowledge about the function. Isaac Newton (1642- 
1727), Joseph Raphson (1648-1715), Leonhard Euler (1707- 
1783), Joseph Lagrange (1736-1813), Carl David Tolme 
Runge (1856-1927), Martin Wilhelm Kutta (1867-1944), 
George Dantzig (1914-2005) and so on have developed the 
numerical differentiation. 

Programming language C is very flexible and powerful. It 
was originally designed in the early 1970s. It allows us to 
maximum control with minimum command. It is recognized 
worldwide and used in a multitude of applications especially 
in Numerical solution of Differential Equations. Along with 
other numerous benefits, we have used programming 
language C in this paper. 

The outline of this paper is as follows: Section 2 contains 
the description of the existing methods of numerical 



differential equations with methodology. In Section 3, we 
develop a numerical simulator, using the programming 
language C, which gives us the solution of an ordinary 
differential equation simultaneously regarding four popular 
existing methods namely Euler, modified Euler, Runge- 
Kutta second and fourth order. Moreover, the technique 
identifies the method that gives the best solution comparing 
with possible exact solution of the problem. We devote 
Section 4 to a trial of the simulator for a specific problem. 
Conclusions are given at the end at Section 5. 

2. Existing methods 

We give a brief description of the existing methods of 
numerical differential equations like Euler, modified Euler, 
Runge-Kutta second and fourth order in this section with 
their methodology. 

2.1 Euler method 

In mathematics and computational science, the Euler 
method is a first-order numerical procedure for solving 
ODEs with a given initial value. It is the most basic explicit 
method for numerical ODEs. 
We consider the differential equation 



[iii] 



y' = f(x,y) (l) 

with the initial condition y(x ) = y (2) 
Suppose that we wish to solve the equation (1) with (2) for 
the value of yat x r = x + rh (r = 1,2,3,...). 

Integrating (1) with respect to X from y to y l and X to 

Xj , we get 



y = y + \f(x, y)dx 



(3) 



Assuming that f(x, y) = f(x ,y ) in X < X < X x , this 
gives Euler's formula 

y { ~ y +hf(x Q ,y ) [for x-x Q =h] (4) 

Similarly for the range X i < X < X 2 , we have 



7i =yi + ]f(x,y)dx 
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Substituting /(x 1 ,y 1 )for f(x,y) where X l < X < X 2 , 

we have y 2 « y t + h f(x lt y^) [since X 2 — X x = h ] . And 
proceeding in this way, we obtain the general formula 
y n+ i *y n +hf(x n ,y n ), n = 0,1,2,... 

2.2 Modified Euler's method 

Integrating (3) by means of trapezoidal rule to obtain 

y, =y +(h/2)[f(x Q ,y ) + f(x l ,y 1 )] (5) 
We thus obtain the iterative formula 

= y + (h/2)[f(x ,y ) + f(x 1 ,yl n) )] ; 
n = 0,1,2,... (6) 
Where y l is the nth approximation to y l . The iterative 
formula (6) can be started by choosing y x from Euler's 
formula y° = y + h f(x , y ) 

2.3 Runge-Kutta method (second order) 

Substitute ^ = y +hf(x ,y ) on the right side of 
equation (5), we obtain 

y 1 = y + (h/2)[f + f(x +h,y +hf )] (7) 

Where f = f(x , y ) andx 1 = X + h . 

Now set k l = h f and k 2 = h f(x + h, y Q +k x ) so 

that equation (7) becomes y x = (1/ 2)[k x + k 2 ] ■ 
This is known as Runge-Kutta second order formula. 

2.4 Runge-Kutta method (fourth order) [iv] 

We mention the fourth order formulae defined by 

y i =y + W A + W 2 k 2 + W 3 k 3 + W 4 k 4 (8) 

Where 

K = hf(x ,y ) 

k 2 =hf(x Q + cc h, y Q + fi k x ) 

k 3 = h f(x + <x x h, y + f} x k x + vjcj 

k 4 = hf(x + a 2 h, y + f3 2 k x + v 2 k 2 + S x k 3 ) (9) 

Where the parameters have to be determined by expanding 
both sides of (8) by Taylor's series and securing agreement 

7 4 

of terms up to and including those containing n . The 
choice of the parameters is, again arbitrary and we have 
therefore several fourth order Runge-kutta formulae. If for 

example we set CC = f3 = CC l = 1 / 2 , CC 2 = 1 , 
/?, =(V2-l)/2, /? 2 =0, v 1= (V2-l)/2, 
v 2 = -l/V2, S l = (V2+l)/2, W y =W A =\l6, 
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W 2 =(V2-1)/3V2 and W 3 =(V2+1)/3V2, We 
obtain the method of Gill, whereas the choice 

a = a 1 =1/2, ^=^=1/2, p x = J3 2 =v 2 =0, 

a 2 =V l =1, W, =W 4 =l/6and W 2 =W 3 =2/6 leads 
to the fourth order Runge-Kutta formulae, whereas 
fcj = hf(x ,y Q ) 

k 2 = h f(x +h/2,y +k l /2) 

k 3 = hf(x + h 1 2, y + k 2 1 2) 

£ 4 = hf(x +h, y + k 3 ) 

Consequently, y x = y Q + (l/6)(k l + 2k 2 + k 3 + k 4 ) 
3. Algorithm 1 " 

We here give the algorithm that can solve differential 
equation numerically regarding four methods- Euler, 
modified Euler, Runge-Kutta second and fourth order 
simultaneously. Moreover, it will ensure which method 
gives the best solution with comparing exact solution. 

INPUT: function f(x, y) , initial condition (x Q , y Q ) , 

interval h , value of X , direct result r . 

Step-1: Set n = (x- x )/h , y 00 = y , y oe = y , 

y 02 =°' y 4 =°' yi = y m + hf(x Q ,y 00 ) 

Step-2: Set i = 1 

Step-3: While i < n repeat Step-4 to step-7 
Step-4: Set j = 1 
Step-5: While j < /repeat step- 10 
Step-6: Set x x = x + h , 

yw = yoo+( 1/2 ) h f( x o>yo)+f( x i>yi)> y\ =yw 

Step-7: Set y u = y 0e + h f(x Q ,y 0e ) , 

Ki= h f( x o>yo4)> 

k 22 = hf(x +h/2,y m +k n /2), 

k 33 = hf(x + hi 2, y Q4 + k 22 1 2) , 

£44 = hf(x Q +h,y 04 +k 33 ) , 

y x r 4 = y 04 +(k n + 2k 22 + 2k 33 +k 44 )/6, 

k x =h f(x ,y 02 ) ,k 2 =h f(x + h, y 02 +k x ), 

yi r 2 = yo2 + ( k l + k 2) /2 > y e = yu> ym = y l0 > 
y 4 = ys 4 , y 02 = ys 2 , x = x n . 

Step-8: Set 

a = \yu- r \' b = \y l0 -r\, c = \y x r 2 -r\, 
d = \y { r 4 -r\ 
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If ( a < b and a<c) 
lf(a<d) 

yu 

Else 

yi r 4 

lf(b<a and b<c) 
lf(b<d) 

y w 

Else 

yi r 4 

If ( c < a and c < b ) 
lf(c<d) 

yi r 2 

Else 

yi r 4 

Output: y le , y w , y { r 2 and y^r 4 with message which method 
gives best solution. 

4 Findings 

We are now offering a trial of the simulator for the 

dy 

numerical differential equation — = — y , with condition 

dx 

y(0)=ltofind y(0.04). 
4.1 Input 

This is a program of differential equation. 
Enter Equation: 

_ dy 

dx 

Enter value of h : 
0.01 

Enter value of X : 


Enter value of y : 
1 

Enter target: 
0.04 

Type the exact value 
0.9607894 



4.2 Output 

The required interpolated value of y[0.0400] in Euler 
method: 

The required interpolated value of y[0.0400] in modified 
Euler method: 1 

The required interpolated value of y[0.0400] in Runge- 
Kutta method: 0.960790 

The required interpolated value of y[0.0400] in Runge- 
Kutta 4 th order method: 0.960789 

The Runge-Kutta 4 th order method is the best. 

5. Conclusion [v] 

In this paper, we develop a simulator incorporated with the 
traditional Euler, modified Euler, Runge-Kutta second and 
fourth order methods for solving ordinary differential 
equations. It been observed in the paper that the result 
obtained according to our procedure is much nearer to the 
exact solution. We therefore, hope that this algorithm can 
solve differential equation problems numerically and can 
save our time and labour. Finally, we have seen that Runge- 
Kutta fourth order gives the best solution among these. 
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